Spectral Functions of Local Operators In The 
Tomonaga-Luttinger Model 



S. P. Strong 

Theoretical Physics, University of Oxford, 1 Kehle Rd., Oxford, 0X1 3NP, UK 

and 

Joseph Henry Laboratories, Princeton University, P.O. Box 708, Princeton, NJ 08544-0708 

(May 20, 1994) 

Abstract 

I give a simple general prescription for computing the spectral functions of 
local operators in the Tomonaga-Luttinger model from the space-time correla- 
tion functions. The method is significantly simpler than directly transforming 
the space-time Greens function and allows a physical interpretation of the sin- 
gularities encountered in the spectral function. 
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Typeset using REVTgX 



The Tomonaga-Luttinger model provides an fascinating example of a soluble model 
of an interacting Fermion problem with a non-Fermi liquid groundstate. The low energy 
eigenstates of the system consist of bosonic degrees of freedom representing collective modes 
of the Fermi surface There are no single particle, Fermionic, low energy eigenex- 

citations. Correspondingly the spectral functions of local operators in the model display 
unusual features including multiple singularities with non-trivial power laws. The proper- 
ties of the single-electron spectral function were studied by Meden and Schonhammer |^ 
and Voit by Fourier transforming the single-electron Greens function. This procedure 
is entirely correct; however, it is somewhat unintuitive and the results take on complicated 
form for the case of models with spin. I give here a prescription for calculating the spec- 
tral functions of local operators from the appropriate space-time Greens functions based 
on the observation that many of the complications encountered in the Fourier transforming 
of the relevant Greens functions arise from the fact that a general local operator in the 
Tomonaga-Luttinger model involves the creation of four distinct types of bosons and these 
complications are circumvented when the independent nature of the different boson types is 
considered. 

Using forms suitable for the calculations of spectral functions in the low energy, long 
wavelength limit, the single electron creation li^)) and annihilation {i/jr or op- 

erators, the local singlet {ipji^i^{x)'i/jL^i{x)) and triplet {'ipR,\{-'^)'^L,i{x)) pairing operators and 
the local charge (-^Jj +■0]^, and spin {iIjI^^^{x)'iIjl,i{x) + 'il!\^^-^{x)ipR^i{x)) 

density wave operators can all expressed as sums of operators of the form : 



ttj) = exp(z ^ aj(t>'j{x)) exp(z ^ aj(t>j{x)) 



(1) 



j j 



where 





where a^ j creates a density wave excitation of the j type fermions of the non-interacting 



model: 
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Y.i^]{k + q)i^,{k) (3) 



The i index runs from one to four with with one corresponding to the left-moving, up-spin, 
two to right-moving, up-spin, three to left-moving down-spin and four to right moving down 
spin fermions. For the non-interacting Luttinger model (and in the low energy limit of the 
Tomonaga model the ^ create eigenexcitations: 



^0 = VpY.i-iy^{{-iyq)qal^a,, (4) 



In the low energy limit of the interacting the Tomonaga-Luttinger model the eigenexcitations 
are linear combinations of the a\ , 

H^nt = Y.{-^yQii-^y^)mb\A^3 (5) 

where 

= E PA a,,, + — 4^4,,) (6) 

k 

The original operators can then be expressed in terms of the new boson operators: 

$(x;7i) ~ exp{iJ2ljCj{x)) exp{iJ2ljQ{x)) (7) 
j j 

where 

k 

SO that 

7j = E PimO^rn (9) 
m 

where /?~^ is the inverse of (]. 

The operator, $, is a product of operators coupling to the different species of bosons, 
resulting of in a product form for the asymptotic, space-time Greens function H,|5|J^,p!0[] : 



(10) 



There is some ambiguity in the choice of /3's required to diagonahze the Hamiltonian. The 
choice: Pi or 2,i = Pi or 2,i±2 and ^3 or 4,i = —Pi or 2,i±2 IS always possible and results in 
Pi = -yf and, for that case, c and s subscripts refer to charge (symmetric combination of up 
and down spin bosons) and spin (anti-symmetric combination of up and down spin bosons) 
bosons, respectively. 

The analytic properties of G are quite complicated and Fourier transforming directly to 
find the spectral function is generically rather involved even numerically 0. However, 
on physical grounds one expects a significant simplification can be achieved if on recognizes 
that the form of G is a simple product resulting from the the fact that $ was a product 
of independent operators. The spectral function of the product of operators of the form 
exp{i J2j IjCji^)) should therefore be given by a convolution of the spectral functions of the 
individual operators: 

1 /"^ r 

p^{k,u) = Re- dte'''' dxe~'''^GS\<^\x,t)^iO,0)\GS) (11) 

TT JO J 

= Y.\{mmO, 0)\GS)\^6{uj - Em)S{k - p^) (12) 

m 

= E (niKlexp(^E74(^))|G''5)pV^^"^^'"^)'^^^-^^'"^) 

mi,m2,m[i,m4 ^ j j J j 

= / n d^jdkj pj (kj ,Uj) S{k-Y.^j)^i^ -J2 ^i) (14) 
j 

where 

Pj{kj,Uj) = E K"^il exp(iE7jCi(2;))|G'>S)|2<5(^j - Em^)6{kj - pm^) (15) 

1 roc f 

= Re- dte'^' dxe-'^''{GS\exp{-iJ2ljQiix,t)exp{iY,ljCj{i^'^)\GS) (16) 
^ ° j j 

The utility of this decomposition for the Tomonaga-Luttinger model results from the 
simple form of the spectral functions of the various parts of $. The spectral function for the 
operator exp{iJ2j IjCji^)) is given by: 

p,{k, u) = \k\-^-P^5{uj - (-l)i^; .A;)e(u;)e-l'=l/'^^^^^^^^r(l + Pj)C{A) (17) 
where C(A) is a cutoff dependent constant and vj is the velocity of the type j bosons. 



The convolutions are thus significantly simplified by the presences of multiple delta func- 
tions. For non-zero 7's there are delta function of this type plus two overall delta func- 
tions for 2N variables of integration. Thus in the worst case the spectral function involves 
two non-trivial integrations, but even these are only of simple power laws. 

For example, for the spinless Tomonaga-Luttinger model the spectral function for the 
electron creation operator involves only one type of right and left moving boson and there 
are no non-trivial integrations to perform. The piece of the spectral function coming from 
the right Fermi point is given by: 

p{k,uj) (X / duJidu2dkidk2S{u — ui — U2)5{k + ki — k2)S{uJi — vki)d{uj2 — vk2)ki ^^k2 (18) 



where k is measured from kp. The singularity at the onset frequency arises because at that 
frequency it is possible to insert very little momentum into the left moving bosons and the 
matrix element for this is singular as k ^ 0. 

This is the general source of all power law singularities which arise in the spectral func- 
tions of local operators in the Tomonaga-Luttinger model: the divergence of matrix elements 
to create states with very little momentum carried by one or more types of bosons. For in- 
stance, including spin with spin-dependent interactions so that the electron creation operator 
coupled to four different kinds of bosons would result in a singularity at = Vck with expo- 
nent — 1 — pi — p3 —pi, resulting from the convolution of the singularities contributed by the 
three species of bosons for which u 7^ Vck. (Here the jo's are those entering the real space 
single electron Green's function and may be obtained from the Bogoliubov transformation 
which diagonalizes Hint once all interactions have been written in bosonized form. [|ll|) The 
piece of the creation operator spectral function coming from the right Fermi point is given 




(19) 



by: 




00 



duJiduj2duj^duJidkidk2dk^dk4^ 5{uji — Vcki)8{uj2 — ^£^2)^(^3 — i^s^3)^(^4 — Vgk^) (20) 



5{uj — uji — U2 — uj^ — Ui)5{k + ki — k2 + k-s — k4)k 
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-Pl J,-l-P2 t,-l-P3 

2 3 4 
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oc / duiduj2Q(uJi)e(uj2)e(uj - Vsk ~ (1 + —)uji - (1 - —M (21) 

JO Vc Vc 

Qiu + vsk - (1 - - (1 + -)u2) 

Vc Vc 

L0^'-P'uj2'~'''{uJ - Vsk - (1 + '^)ui - (1 - ^)a;2)-^+^^ 
+ y^k - (1 - ^)uJi - (1 + ^)a;2)"'+P* 

Vc Vc 

which has exactly the singularity at u; = Vck expected. In general a singularity in the above 
can occur for u = Vik and the exponent of the singularity is always given by —1 — J2j=iiPj- 
Note that Eq. |2l] gives the low energy, long wavelength spectral function for any local 
operator whose real space Green's function is given by Eq. |10|, so the above formula gives 
the spectral function for any local operator in terms of the exponents and velociites entering 
into its space-time Green's function. The same type of singularity at ci; = Vik is present for 
all operators with the same sum of pj's for j ^ i. It is only through the dependence of this 
and other such sums on the actual operator that different operators acquire different types 
of singularities. 

This method gives a simple, intuitive way of obtaining the spectral function of all local 
operators in the Tomonaga-Luttinger model, however it suffers from two drawbacks. The 
first drawback to this approach to obtaining the spectral function arises from the fact that 
the inclusion of the cutoff in Eq. |10| is not strictly correct. The spectral functions obtained 
from G will therefore not correctly reproduce cut-off related features and result, in 

general, not properly normalized. Once properly normalized they should properly reproduce 
all of the low energy physics. 

The second drawback is that the real part of G{k,uj) is not obtained directly in this 
approach, but requires a further Hilbert transform of the spectral function. Since all of the 
physical information is encoded in p{k,uj), this is not that serious a problem, but for nu- 
merical applications where the real part of G is required it will be a significant disadvantage 
for this approach compared with methods based on direct Fourier transformation of the real 
space Green's function. 
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